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Abstract 

In this paper, I consider a thin suspended liquid film, surrounded by 
a different fluid. Examples of such a system are soap films and liquid 
crystal films, surrounded by air. They are considered good models for 
two dimensional fiuid dynamics, although air drag is known to introduce 
strong deviations. I exactly solve the hydrodynamic equations of the sys- 
tem composed by the liquid layer and the surrounding fluid in order to 
evaluate the correlation function of the thermally excited velocity fluctu- 
ations. Both the temporal and the spatial statistical properties are the 
ones of a two dimensional, hypotetical fluid, only for temporal and spa- 
tial frequencies higher than critical values; at lower frequencies they show 
signiflcative deviations due to the interactions of the liquid layer with the 
surrounding fluid. 

1 Introduction 

A thin, freely suspended liquid layer is generally considered a system that fol- 
lows two dimensional hydrodynamics, on length scales longer than the layer 
thickness. Many experiments have been performed on soap films, that consist 
of two monolayers of amphiphilic molecules around a layer of water. The layer 
thickness ranges from the 4nm of the Newton black film to many micrometers; 
features down to a fraction of a millimeter should be considered two dimensional 

Also smectic A liquid crystals can form films, by a different effect. It spon- 
taneously form a layered structure, with the elongated molecules alligned to 
the surface normal. When a film is drawn from this kind of liquid crystal, it 
is formed by a set of monolayers, ranging from two to many hundreds [T]; the 
thickness can be selected with extreme accuracy. Because of the absence of 
translational order in the monolayers, the film behaves as a two dimensional 
system. 

Many papers [6l [71 [TOj HI [5] describe experiments in which turbulent flows 
have been detected, and discuss the agreemen with two dimensional hydrody- 
namics. One result is that the role of the air surrounding the film cannot be 
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neglected. In turbulent flows, it damps the turbulent vortices. From a theoret- 
ical point of wiev, an exact evaluation of this effect is extremely difficult; for 
example, Rivera et al. modeled this damping as a linear drag term in the two 
dimensional Navier-Stokes equations, in order to describe velocity fluctuations 
observed in a turbulent flow. They evaluated the amplitude of the drag term 
from experimental data; from this value, they found the dissipation due to air 
friction to be a significant energy dissipation mechanism in their experimental 
system. 

Another well known effect of air friction on thin liquid films can be found 
in diffusion of particles, whose diameter is about the thickness of the layer. 
Neglecting air friction, the mobility of a particle in a two dimensional system 
should be infinite, that is, a particle subjected to a steady force should accelerate 
indefinitely. This would result in an infinite diffusion coefficient. Saffman [llj 
found that air friction reduces the mobility to a finite value: he evaluated the 
diffusion coefficient and found that it strongly depends on the viscosity of the 
surrounding fiuid, and diverges as the viscosity vanishes. Saffman formula has 
been tested both in soap films [3] and in liquid crystal films [1 . 

The diffusion coefficient can be evaluated once we know the correlation func- 
tion of the thermally excited velocity fluctuations; this means that the velocity 
correlation function must be strongly affected by air friction. In this paper, I 
evaluate the correlation function of the velocity field of a thin suspended liquid 
film, by solving the hydrodynamic equation of the whole system, including the 
liquid film and the fluid surrounding it. 

Then I discuss the result, showing that both the temporal and the spatial 
statistical properties are the ones of a two dimensional, hypotetical fluid, only for 
temporal and spatial frequencies higher than critical values; at lower frequencies 
they show signiflcative deviations due to the interactions of the liquid layer with 
the surrounding fluid. 

Last, I derive approximately the value of diffusion coefficient, and I compare 
it with Saffman formula. 

2 Derivation of the power spectrum of the ve- 
locity fluctuations 

I describe the film and the fiuid surrounding it by using Landau's fiuctuating 
hydrodynamics equations: 

/ Pom'^o {x, t) = r^oV^-uo {x; t) - Vpo {x, t) + ^ [^lii {x, z = 0,t)- (x, z = 0, t)] 

\ Pl-§^Ul,2 {X, Z,t) V^Ml,2 {X, Z,t)~ Vpi,2 (X, Z, t) 

(1) 

where ui {x, z, t) and U2 {x, z, t) are tree dimensional vector fields, describing the 
velocities in the space above and below the film; uq {x, t) is a two dimensional 
vector field, describing the velocity in the film; 770 and ryi are the shear viscosity 
of the film and the surrounding fluid; po and pi are the densities of the fllm and 
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the surrounding fluid, h is the thickness of the film; pi 2 {x, z, t) and po {x, t) are 
the pressures. 

I suppose that both the fluids are incompressible, the film is placed dX z — Q, 
and the fluid does not slip at the surface of the sheet: 

V-uo(.^,t) = (2) 
V ■ ui^2{x,z,t) = (3) 
ui^2{x,z^Q,t) = [uo(f,t),0] (4) 

The pressures pi^2 [x, z, t) and po {x, t) must be determined in order to fulfill the 
incompressibility equations Eqs ([U [3]) . 
In Fourier space: 

{-ipoUJ + rjaq^) uq [q, uj) = iqpo {q, uj) - J [ui {q, qz,uj)- U2 (<f, g^, w)] 
[-ipiw + 7/1 + g2^] ui^2 {q,qz,uj) = i [q,qz]pi,2 {q,qz,(^) 

(5) 

with the conditions: 



q-UQ{q,uj) = (6) 

[q,qz]-ui,2{q,qz,^) = o (7) 

ui.2 {q,qz,i^)dqz = [wo(g,a;),0] (8) 



where the Fourier transform is defined as follows: 

•^('^^ = 7:7^ / ^"^"^fi^)^^ (9) 
\2tt) j 

Now I use the incompressibility equations Eqs ((21 El in order to evaluate the 
pressures: 

f (-ipow + 7?of?^)'"o(g,w) = -«ir (^^9^) S q^[^i{q^qz,^) - U2{q-,qz-,u})]dqz 
\ [-ipiu + 771 (g2 + qiy\ Ui^2 {q, qz, = 

(10) 

with the conditions: 

q-Uo{q,LL>) = (11) 
K^z] • Wi,2 (<f,gz,w) = (12) 

ui,2{q,qz,^)dqz = [uo{q,uj),0] (13) 
Then, I express the vectors on a suitable basis: 

\ 

\Jql + ql 



no {q,Lo) = uoK(g,^)4^i+^OT(g,^)%=ji (14) 



ql + ql 
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Wl,2(g,gz,w) = Ul,2R{q,qz,UJ) +m,2T(g,gz,t^) ■ " 



-zpiw + 7?i ( + g2 



- i/Oio; + rji [q^ + 



\Jql + ql + ql yix + ql 

1^ X [-qxqz,-qyqz,ql + ql] 
m,2z(g,g^,t^) , , = (16) 

^Jql + ql^Jql + ql + q^ 

The subscripts R and T label the radial and tangential components of both 
and Ui,2, and Z labels the component of wi,2 perpendicular to the film. The 
equations become: 

{_-ipoUJ + rioq^)uQT{q, w) = -ix / 1^ ["i^ (^^ 9^. - '"2T (g, g^, w)] dq^ 
Wi,2T {q,qz,i^) = 
ui,2Z {q,qz,oj) = 

WoR (g, w) = 
Ui,2_R {q,qz,i^) = 
/ Ui,2T {q,qz,uj)dqz = uoT (g,w) 
/ Ui,2Z (g, g^, w) Jl^dq^ = 
. /ui,2z(g,g^,w) |jJ|^rfg^ = 

(17) 

Only the components uqt and iii,2T are needed in order to evaluate uq; in 
the following, I will drop the subscript T. 

(-ipouj + rjoq"^) uq {q,uj) = ^ q^ [m {q,qz,u}) - U2 {q,qz,uj)]dq^ 

[-ipiu + r^i (g^ + g^)] ?ii,2(g,gz,w) = (18) 

{q,qz,oo)dqz = Uq (g,w) 
In real space: 

^uo (g*, t) = -i^oq^uo {q, t) + ^ (q, z = 0,t) - (g, z = 0, t)] 

^ui,2 (g, z,t) = -i^ig^Mi,2 (g, z,t) + ui-^ui^2 (g, z, t) 

Ul,2 {q-,z = 0,t) = Uq {q,t) 

(19) 

where v = r]/p. 

In order to simplify the following calculations, I impose periodic boundary 
conditions on z, with a period L; later, I will evaluate the limit for L — >■ oo. I 
express the hydrodynamic variables on a complete orthonormal basis: 



u,,2{q,z,t) = Y.{<it)fni^) + bi{t)f^iz)}+cHt)f{z) (20) 

n=0 

oo 

= Y.<it)fniP) + cHt)fiP) (21) 
n=0 

The eigenfunctions are: 

f^{z) = N^s\n[A^{z-L/2)] (22) 
f^{z) = N^coa[B„{z-L/2)] (23) 
f(^{z) = N'^cosh[C{z-L/2)] (24) 
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where the wavevectors, for n > 0, are: 



Bn is implicitly defined by: 



2 pi Bn 

and the only value of C is implicitly defined by: 

-(2 



2>i C 
The normalization constants are: 



l^Pol =tanhfc^ 



= 1/1 



2/ipo 



1/^ 



1 



^ 2/ipo 



cos^ ( -B. 



= 1/^ 



Lpi 



\ 777^ + 77 cosh^ C 



The eigenvalues are: 



'■n 



a:: = -ui (g2 + A: 



-u, {q^ + El) 
A^ = -^i(<?^-C2) 

The eigenfunctions are ortogonal with respect to the scalar product: 

Pi 



f-9 = 



hpo 



f{z)g{z)dz + f{0)g{0) 



The evolution equations for the coefficients are: 

f (t) = XM (t) 
MW = Af6f(t) 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 
(32) 
(33) 



(34) 



(35) 



In order to use equipartition theorem, I evaluate the kinetic energy of the 
fiuid: 

E=^Poh j \uo{x,t)\^ dx +^pi j \ui{x,z,t)f dxdz (36) 
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In Fourier space, the energy can be decomposed as the sum Et + Ez^ where: 



Et — 



-^—pQti / \uQ{q,t)\ dq-\ — - 



Pi J \ui{q,z,t)f dqdz (37) 



The integrals can be expressed in terms of the scalar product defined in Eq. 



Ex 



poh / u {q, z,t) ■ u {q, z,t) dq 



Due to the orthonormality of the eigenfunctions I used: 



(27r) 



[n=l 



\bff 



dq 



From equipartition theorem: 



5n,7nS {q + cf) ^^'^f 



(27r) Pah 
KbT 



(27r) Poh 

From Eq. (1351) I can derive the time spectrum: 



c«Hc« {uj') 



d (q + q ) (LU + LJ ) 



(38) 
(39) 

(40) 
(41) 

(42) 
(43) 



Using Eq. (j2ip . I evaluate the power dynamic spectrum of the velocity fluctua- 
tions of uq: 

+ 00 ,g 



—A,, 



n=l 

KbT -X^ 



\B2 



iV-'^cosh'' I C- 



L 



where the dynamic power spectrum is defined as follows: 

(uo (<f, 1^) Uq ((f, w')) ^S{q + q')S{uj + u') Su {q, uj) 



(44) 
(45) 

(46) 



Now I use Eq. (|26p in order to express any trigonometric function of i?„ and 



C: 



KbT g2 ,^i{q' + Bf, 



KbT i^i(q^-C^) 



2 
C2 



5„(g,w)= (47) 
+ (48) 



An^poh oj^ + „f {q^ - C^f 



T J _£l /^2 1 /iPo 



C2 



C2+ 1 



(49) 



^1 . 
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From Eq. 

Bn+l — B„ 



27r 



c = 




1-^ 

Vl 



Pi 



(50) 

(51) 



For L ^ oo: 



KbT 1 



+ 00 



!/i (g2 + 



vl (g2 



hpo 

KbT 



1 hpo 
4 pi 

i^l (g2 - 



5„(g,c^) = (52) 
-dB+ (53) 



1 



C72 



(54) 



I define: 

a 
90 

h {x, y) 

l{x,y) 
A{x) 
and obtain: 
Su (g, w 



1 

Pi 

4 

TT 



1^0 

+ 00 



1 



B' 



B' 



y 



B^y 4B2 + (Q,a;2 ^ 52 



1 — a V 

r2(x) 



1 



ny"^ + (x) 
1 



1 



ax^ 



KbT 



1 



2M^, 



(55) 
(56) 

(57) 

(58) 

(59) 

(60) 

(61) 



The term with I (x, y) has a Lorentzian behaviour, and represents the sponta- 
neous velocity fluctuations of the film, damped by the surrounding fluid. The 
other one describes the velocity fluctuations induced by the fluid on the fllm. 
It's worth noting that the only spatial lengthscale is 1/qo, while two different 
temporal lengthscales are involved in the two processes: i^oQa for the fluctuations 
of the film, and viQq for the velocity fluctuations of the surrounding fluid. 
Some explicit results can be obtained by evaluating the integral: 



+00 



1 



52 



1^2 452 _^ (2;2 



B^y 



rdB 



1 



xVjF+T 2^\Jx^ + 1 + (a + 1) X 

(62) 
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With a = 0, I obtain, as expected from the equipartition theorem: 



+ 00 j^^rj. 

Su {q, duj = — — - (63) 

The power spectrum on a; = can be evaluated by means of Eq. (1621) : 

S^{q,u; = 0) = \ (64) 

i-K-^rjoh q (qc + q) 

where: 

9C = 7— , (65) 

Kb is the Boltzmann constant, T the temperature, etag and etai are the shear 
viscosity of the film and the surrounding fluid, and h is the thickness of the film. 

The power spectrum has a q~^ dependence for q ^ qc', at smaller wavevec- 
tors, the divergence saturates to q^^, due to air damping. For a soap film in 
air, the values of densities are po = lO'^Kg/m^ and pi = l.SKg/m"^. The ef- 
fective viscosity of the film is I'o = 1-6 • fO^^m^/s, [9], evaluated by means of 
Trapeznikov relation [T^]. The viscosity of air is h'l — 1.43 • 10~^m^/s [S]. From 
these values, for a 2pm thick film, qc ~ 1.2 • lO'^m^^, corresponding to a 0.5mm 
wavelength. 

From the power spectrum of the velocity fluctuations, we can derive the 
diffusion coe5icien[2]: 

D^n [ Su{q,uj^O)d\ (66) 
Ja 

where the cut off wavevector A is of the order of the inverse of a, the radius of 
the diffusing particle. By integrating the static power spectrum of Eq. (|64|) : 



KgT 

D= " 



2'KrjQh 



In In (2) 



(67) 



This is the well known Saffman formula for the diffusion coefficient of a particle 
on a liquid film [TTl El [1] . 

This result can be compared with a phenomenological theory [9]. 



3 Time correlation function 



By Fourier anti-transforming Eq. ()6ip in cj, I obtain the time velocity time cor- 
rlation function: 



Su (q, t) 



KbT 
An^poh 



where: 



1^1 



(68) 



(69) 
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= - / e-y^ ^dB (70) 

Trio 4B2 + (^2+52)2 

= vo + {v^-vo)2^^^±^^ (71) 



go 



(72) 



V 1 + a;^ 

<Zo = (74) 

«Po va 
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